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Abstract In this paper we give a general introduction to quantum critical phenom- 
ena, which we practically illustrate by a detailed study of the low energy properties 
of the spin boson model (SBM), describing the dynamics of a spin 1/2 impurity 
(or more generically a two-level system) coupled to a bath of independent harmonic 
oscillators. We show that the behavior of the model is very sensitive to the bath spec- 
trum, in particular how the properties of the quantum critical point in the SBM are 
affected by the functional form of the bath Density of States (DoS). To this effect, 
we review the renormalization group (RG) treatment of the SBM for various bath 
DoS, based on an unconventional Majorana representation of the spin 1/2 degree 
of freedom. We also discuss the derivation of Shiba's relation for the sub-ohmic 
SBM, and explicitely derive an effective action vindicating the quantum to classical 
mapping. 



1 Introduction 

Quantum Phase Transitions (QPT) have recently become a widespread topic in the 
realm of modern condensed matter physics. QPT are phase transformations that oc- 
cur at the absolute zero of temperature and are triggered by varying a temperature 
independent control parameter like pressure, doping concentration or magnetic field. 
There are various examples of systems showing quantum critical behavior, which in- 
clude the anti-ferromagnetic transition in heavy fermion material like CeCu6_ x Au x , 
that is brought about by changing the Au doping ifTTl . Another prototypical exam- 
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pie of a system exhibiting quantum critical behavior is the Quantum Hall Effect, 
wherein a two-dimensional electron gas is tuned, via an externally applied magnetic 
field, through a quantum critical point (QCP) that intervenes between two quantized 
Hall plateaux. Other examples of QPT include the ferromagnetic transition in metal- 
lic magnets as a function of applied pressure, and the superconducting transition in 
thin films. 

Since there are such a wide range of experimentally accessible systems that show 
quantum critical behavior, it is imperative that we understand QPT at a fundamental 
level. We shall here endeavor to do just so by giving an introductory account of this 
fascinating phenomenon. As a striking illustration, we will be comparing and con- 
trasting QPT with the case of more usual thermal (classical) phase transitions (as 
will be seen later on, thermal phase transitions are also referred to as classical tran- 
sition, since quantum fluctuations become unimportant in their vicinity). Let us first 
begin by discussing the ferromagnetic transition, in order to better illustrate the rich 
phenomenology of phase transitions (both classical and quantum). The route that 
we take here to understand the fundamentals of QPTs is as follows: We shall first 
review the basic phenomenology of classical (thermal) phase transitions. Then, we 
shall illustrate via heuristic arguments how quantum fluctuations can be disregarded 
in the vicinity of a thermal phase transition. These arguments also provide clues to 
the domain in the phase diagram where one expects quantum fluctuations to dom- 
inate. Also, we shall briefly discuss the question of observability of QPTs. Finally, 
we shall end with a discussion of the so called quantum to classical mapping. 

Let us first start with classical (thermal) phase transitions. As a physical system, 
say a ferromagnet, approaches its ordering, there is a length scale called the corre- 
lation length, 4 that diverges in a power-law fashion when one comes closer to the 
critical point, £, ~ \ t |~ v , so that the system becomes progressively self-similar. Here, 
t is dimensionless parameter characterizing the distance to criticality, t = — ^ (for a 
thermal transition), and T c is the critical temperature where the phase transformation 
occurs. Now, the above divergence of the correlation length encapsulates the infor- 
mation that the fluctuations of the order-parameter (say the magnetization) become 
spatially long-ranged as the system approaches the critical point. Analogous to £ 
one can define a time scale § T , that also diverges as a power law as one approaches 
a second order transition. Thus, we have: 

k=s<=i*r vi a) 

The quantity z that controls the divergence of <i; T is the so called dynamical exponent. 
Now, associated with this time scale we can define a frequency scale co c 1 /^ T , and 
through it a corresponding energy h(0 c , which encodes information pertaining to the 
energy scale related with order-parameter fluctuations. This quantity hco c competes 
with kaTc, the typical energy associated to thermal fluctuations. Now, the question 
of importance of quantum fluctuations can be re-cast into a query of which among 
these two energy scales prevails. Since co c — > as one approaches the critical point, 
the energy scale of thermal fluctuations (for any non-zero T c ) always dominates over 
the scale hco c . In other words for a transition that happens at a finite temperature, 
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h(o c <C kgT L . Thus, it can be argued that asymptotically close to a finite temperature 
transition, it is the thermal fluctuations that are the driving mechanism. This irrele- 
vance of quantum fluctuations near a thermal phase transition is the reason to why 
they are given the moniker "classical phase transition". 

Now, from our discussion in the previous paragraph it is but obvious that if the 
transition were to occur at T = (tuned by a non-thermal parameter like doping or 
pressure), then the fluctuations that will drive the transition will be wholly quantum 
mechanical in origin. It is obvious that one then needs to apply ideas from quantum 
statistical mechanics to understand QPT, as pionneered by Hertz in a seminal pa- 
per Q to tackle the problem of quantum criticality in itinerant magnetic systems. 
By using this case of the quantum magnet as a test-bed example, Hertz J7] showed 
that any generic d dimensional quantum system can be mapped onto an equivalent 
d + z dimensional classical model. This statement is referred to as the quantum to 
classical mapping and is of fundamental importance in the field of QPTs. By us- 
ing the quantum to classical mapping one can show that the critical behavior of the 
quantum model is equivalent to that of a classical model but in z higher dimensions. 
Although this mapping is believed to be robust for insulating magnets, it was how- 
ever later shown (T| that Hertz's conclusions were erroneous for a large class of 
itinerant QPT. This break-down in fact occurs due to the presence of soft modes 
in the systems (e.g. the particle-hole excitations in itinerant magnets) other than 
the order-parameter modes. The presence of these modes induces an effective long- 
ranged interaction between the order-parameter modes, thereby altering the critical 
behavior (TJ, as compared to Hertz's original results. 

Since QPT occur at zero temperature, it was initially thought that the study of 
these phase transitions was a mere academic exercise. However, it was soon realized 
that the presence of a zero-temperature critical point (practically inaccessible) can 
actually influence the behavior of the system at finite temperatures. In other words, 
at any finite temperature, the critical singularities associated with the QCP are cut by 
the temperature, so that one observes non-trivial temperature dependence of various 
observables in a so-called quantum critical regime. The calculation of the quantum 
critical regime for various models is well beyond the scope of this work. However 
the interested reader is directed towards the following papers investigating the effect 
of non-zero temperatures on QPT in magnetic systems Il20ll23l . Also, one will refer 
to Sec. I2.2l for a brief description of the quantum critical regime in the spin boson 
model (SBM), the specific model of interest in this manuscript. 

Thus, from the discussion of previous paragraphs, it is obvious that QPT are an 
extremely interesting physical phenomenon to study. As alluded to before, in this 
manuscript we choose to study a specific toy model example, namely the QPT en- 
countered in the spin boson model (SBM), a variant of Caldeira-Leggett type models 
fl9l , wherein a quantum particle is subjected to an external dissipative environment. 
In the case of the SBM, this quantum "particle" is essentially a two-level system, 
such as a spin 1 /2 impurity. While there have been many studies on dissipative quan- 
tum models that focused on the effect of decoherence on intermediate time scales, 
their behavior in the long time limit in the presence of quantum critical points remain 
relatively un-explored. However, the study of such regimes is extremely important 
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as anomalous low energy properties emerge due to quantum critical modes. In other 
words, due to the presence of a QCP, the SBM can display non-trivial dynamics at 
very long times. 

As a more general remark, we note that this model can also be used to study 
quantum criticality at the level of a single spin j impurity embedded in a corre- 
lated system such as Mott insulators ||3][l3][T8), or magnetic metals Rl flOl ). It also 
appears as an effective theory for bulk materials themselves (e.g. in quantum spin 
glasses |fl6ll , heavy fermion compounds |fl5l[T7l ), via the framework of DMFT. 

The SBM is introduced in Sec. [2] The various phases of the SBM and the pos- 
sibility of a QPT between them is discussed in Sec. 12.21 In Sec. [3] we re-write the 
SBM by using a Majorana fermion representation for the impurity spins. Sec.|4]is 
devoted to the derivation of the RG equations by using the Majorana representation 
presented in Sec. [3] In Sec. [5] we look at the consequences of the flow equation 
derived in Sec. [3] Sec.|6]is dedicated to the quantum/classical mapping of the SBM 
to the long-ranged Ising model. Sec.|7]is concerned with the development of a spe- 
cial identity in the SBM model that is used in Sec.|8]to derive the so-called Shiba's 
relation in the case of the sub-ohmic spin boson model. Sec.|9]is dedicated to a dis- 
cussion on the status of the quantum to classical mapping in the SBM, that we use as 
a conclusion and future outlook regarding quantum phase transitions in dissipative 
models. 



2 The Spin Boson Model 




Fig. 1 Fig. [JJis a pictorial sketch of the SBM. It represents a two level system (e.g. a particle in 
two potential minima ±<po, or a spin 1/2 impurity, coupled to a bath of harmonic oscillators (such 
as phonons, nuclear spins,...) via a coupling constant X (wavy line in this figure). 
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As stated earlier, the SBM describes the effect of an external dissipative environ- 
ment on the quantum mechanical evolution of a two-level system. We will introduce 
in Sec. 12. li the general properties of the SBM, and present in Sec. 12.21 its possible 
phase diagram, obtained on heuristic grounds via an analysis of the various limiting 
cases. 



2.1 The Model 

The SBM involves a single spin A impurity S, interacting with a set of bosonic bath 
variables, a,, and a] (in second quantization). The interaction between the bath's 
oscillator displacement and the spin is controlled via a coupling constant A. Thus, 
the SBM Hamiltonian has the general functional form: 

H = -AS X + £S Z + AS z £(a? + «,-) + Y,®i a \ a i- ( 2 ) 

i i 

Here, in Eq.|2] A and £ are the transverse and longitudinal magnetic fields respec- 
tively, applied to the quantum spin. A physical sketch for such a SBM, wherein a 
two-level impurity (when the bias field e is set to zero) is connected to an external 
environment, is depicted in Fig.Q] All that remains to completely specify the model 
is to endow the bosonic degrees of freedom with a spectrum. This bosonic density 
of states (DoS) is taken here to be continuous and power-law like, and conforms to 
the functional form: 

p(co)^8(co-co i )= ( d(co)d(A-co). (3) 

Here, in Eq.|3] A is a high-energy cutoff. When the exponent s is such that < s < 1, 
then the model is said to be in the sub-ohmic regime, the case s = 1 is referred to 
as ohmic, while the case s > 1 is called super-ohmic. In fact, as will be seen in the 
course of this paper, the quantum critical behavior of the SBM is crucially dependent 
on the exponent s controlling the behavior of the bath spectrum. 

It is also convenient to define a new bosonic variable corresponding to the "local" 
displacement: 

<M !>+«•), (4) 
i 

which has an associated DoS given by, 

p^(o) = - is+ A ] ) } al \ gn( ( o)9(A 2 - ( o 2 ). (5) 

To make comparision with existing literature, one can alternatively characterize the 
bath by means of a spectral function: 
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J(co) = ^ttA 2 5(co - O;) = 2na(o s A l - s d((o)e{A - co). (6) 
Here a is the non-dimensional dissipation strength defined as a = (s + l)4j. 



2.2 T/ie P/iases o/f/ie SflM 



Non-universal High Temperature 



Quantum 
Critical 



De-localized 



Localized 



a 



a 



Fig. 2 Fig.[2]shows the generic phase diagram of the SBM model. Here, a quantum critical point 
a c separates localized and delocalized phases, from which a quantum critical region emerges at fi- 
nite temperature. The scale T* is the cross-over temperature below which various physical observ- 
ables revert from quantum critical behavior to those associated with the localized or the delocalized 
phase. 



As promised in Sec. [2] we will study here the possible phase diagram of the SBM 
by looking at situations where either one of the two parameters A or A dominates. 
For instance, let us first consider the case where the dissipative coupling A is set to 
zero. The SBM then becomes equivalent to the case of an isolated spin in a trans- 
verse magnetic field. It is well known that such a system displays Rabi oscillations. 
That is, if one were to start with an initial state pointing "up" along the z direction, 
then the transverse field A periodically drives the system between up and down con- 
figurations. This limiting case A = is in fact adiabatically related to a whole phase 
at non-zero A , dubbed for obvious reasons the delocalized phase, where coherent 
spin oscillations are expected to occur (at least for small enough A). We note that 
the average (S x ) is always non-zero as long as the transverse field A is finite, and 
thus cannot play the role of an order parameter. However, we can pursue a magnetic 
analogy by noting that the longitudinal spin average (S z ) is identically zero in the 
delocalized phase, so that we can really relate this portion of the phase diagram to a 
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ground state with zero magnetization. In fact in the alternative regime, i.e. when the 
dissipation A dominates over the transverse field A, the ground state becomes dou- 
bly degenerate, as can be checked on the trivial limiting case A = 0, where the z spin 
component is clearly conserved within the Hamiltonian (01. A simple physical pic- 
ture emerges, with the system localizing in one of the two minima at ±0o> see Fig-Hi 
Assuming adiabaticity by switching on the transverse field, we arrive to the so-called 
localized phase, wherein the spontaneous magnetization ~ eM^O. Now, 
so far by using heuristic arguments, we have shown that the SBM allows for the ex- 
istence of a de-localized phase, with (<j)\ = 0, and a localized phase, with (<j>) ^ 0. 
Thus, it is quite plausible that a second-order phase transition takes place between 
the two phases. As we shall see explicitely in Sec. [5] there is indeed a second order 
quantum localization/delocalization transition for all < s < 1 . 

The generic phase diagram for the SBM is shown in Fig. [2] where a T = phase 
transition, separating localized and delocalized phases, takes place at a critical value 
a c of the adimensional dissipation strength. The interesting quantum critical regime 
emerges above the critical point at finite temperature, where anomalous behavior 
of all physical quantities is expected. For instance, the longitudinal spin suscepti- 
bility in the quantum critical regime obeys the behavior % Z {T) ~ 1 /T s , as opposed 
to the conventional 1/T Curie-law expected for the whole localized phase. This 
anomalous power-law behavior is a direct signature of the QCP at a c , and will be 
demonstrated in the following sections. We note that for the ohmic s = 1 case, the 
conventional treatment for studying the QPT is to map the SBM into an anisotropic 
Kondo model (AKM), and use previous knowledge on the scaling properties of this 
well-known Hamiltonian. However, in this paper, we shall follow a less well-trodden 
path, namely, performing a renormalization group (RG) calculation directly within 
the SBM, using a spin represention in terms of Majorana fermions (see Sec.[3]for 
further details). This formalism has the advantage that it can be easily adapted to 
perform calculations in the sub-ohmic limit, i.e. (0 < s < 1), see Sec. [4] 



3 The SBM using the Majorana representation 

In Sec. @] we aim to derive the RG equations for the SBM, based on a perturba- 
tive analysis around the localized limit, i.e. A = 0. A technical difficulty on this 
path comes from the fact that the quantum spin i impurity does not follow either 
bosonic or fermionic commutation relations. Thus, standard calculations based on 
Wick's theorem cannot be invoked. One of the many ways to avoid this problem is 
to map the spin--j operator onto fermionic degrees of freedom, which can be done 
in particular using so-called Majorana fermions. The use of this mapping to con- 
densed matter is relatively recent and for a more detailed explaination the readers 
are referred to the following references, Refs. I22ll2p . 

The mapping between the spin A impurity and the Majorana fermions obey the 
following correspondence principle: 
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f = ~t*t (7) 

Here, in Eq. [7] the 77 fields represent a triplet of Majorana (real) fermions 77 = 
(171,172,173). that satisfy the following anticommutation relations {77,-, ?7 ; -} = 8ij. 
Now, in addition to these Majoranas defined above, one can construct another 
fermionic field, <P = 2/771772773, that commutes with the Hamiltonian, and consti- 
tutes hence a conserved quantity, with the constraint <P 2 = j . A very useful relation 
for describing the spin dynamics is given by the correspondence (see ED ): 

if = 20 5^ (8) 

Now, in terms of the Majorana fermions (see Eq. |7), and after the redefinitions 
{S x — > ^3 , S y — > S2 , S z — > Si } which amounts to a n/2 rotation around the y direction, 
the Hamiltonian of the SBM can be expressed as 

H= -i— (771772 -772771) +// B -;A0 772773 (9) 

Now in what follows, we will use Eq. [9] to perform the perturbative RG analysis. 
Before we go on to do so, a word of caution regarding the fermionic mapping is in 
order: Any mapping of the spin j impurity to fermionic operators tends to enlarge 
the dimensionality of the Hilbert space. How such an enlargement is obviated in 
the case of the Majorana representation is technical matter that goes well beyond 
the scope of this manuscript, and the reader is directed to Refs. 12T1I221 for further 
details. 



4 Perturbative renormalization group in the localized regime 



----X----- AA \\ 

Fig. 3 Pictorial representation of the A vertex. 

In this section, we will perform a perturbative RG treatment starting from the lo- 
calized phase, i.e. from the limit in which A = 0. Our plan of action to derive the RG 
equations is as follows: We initially start with a model of free spin (A = 0, a = 0), 
and then perform a perturbative analysis in both A and a, leading to renormaliza- 
tions of the dissipation a and the transverse field A, which depend explicitely on a 
generic cut-off scale A. Following the philosophy of the RG, one aim to compute 
the renormalized parameters at a lower cut-off scale, A', leading to so-called flow 
equations. 



Quantum phase transition in the spin boson model 



9 



The key ingredient in developing the perturbation theory are the free Majorana 
fermion propagator G^ ee , as well as the A vertex shown in Fig. [3] and the A ver- 
tex (first diagram appearing in Fig. The free fermion propagator in Matsubara 
frequency co„ = (2n + \)kT at finite temperature T reads G^ ee (i(Q„) = l/i(O n . One 
can then first construct the vertex function r a related to the dissipative coupling A, 
shown in Fig. [4] The functional form of the vertex function can easily be deduced 
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Fig. 4 Lowest order diagrams involved in the renomialization of the dissipative term X . 



to be 



r a {co 1 A)='- T +'-T^G^{co) 



(10) 



A A 

Here, once again A is cut-off scale, set e.g. by temperature or the bandwidth of 
the bosonic modes, and © is a frequency. The above equation can be effectively 
re-written in terms of an adimensional transverse field h — 2A/A, so that the renor- 
malized dissipation reads: 



r(o),A) = - 



1 exp I 2 In — 

4 F \ (0 



(11) 



Now, as stated in the introductory part of this section, we re-scale the cut-off A to 
A' = A — dA. Under such a re-scaling the vertex function can be re-written as: 



w A ,. A / dA 

r(co,A') = - i + x 



h 2 ( A\ h 2 dA 

1 exp 2 In — H 

4 F \ co 2 A 



(12) 



The above equation can be re-written in the form of the usual RG /3 -function by 
including the frequency dependent vertex function into the redefinition of the cou- 
pling constant. Then in terms of the logarithmic differential d£ = —jr, the Eq.[T2l 
can be re-cast into the form: 



(13) 



dX _ X f 

~dJ ~ ~2* 

and finally more compactly expressed in terms of the dimensionless dissipation a = 

2X 2 /A 2 as 



da 
~d( 



— — = -ah 1 



(14) 
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Fig. 5 Diagrams that are involved in the renormalization of the transverse field A . 



Now, in a similar vein one can calculate the first corrections to the transverse field 
A, with the diagrams depicted in Fig. [5] The technical details of this calculation are 
very similar to the above calculation, and the final flow equation, written in terms of 
the scaled magnetic field h — A/A, reads: 

§ = (!-«)*• (15) 

The RG equations that we have so far derived are for the case of the ohmic damp- 
ing. The derivation of the flow equations in the non-ohmic limit is quite straight 
forward and can be performed by following the technical details elucidated above. 
Thus, for the sake of brevity we will not perform these computations here. Instead, 
we will just quote the results of such an exercise. In the presence of non-ohmic 
dissipation the a flow equation of Eq.[l4]gets modified into 

da , 

— = -ah 2 + (l-s)a. (16) 
at 

However, the flow of the magnetic field h retains its functional form given in Eq.[T5l 
even in the presence of non-ohmic dissipation. 



5 Analyzing the RG flow 



In this section we shall discuss the RG flow equations that were derived above. In 
Sec. 15. II we shall first analyze the j3 functions for the ohmic case (s = 1). Then, in 
Sec. 15.21 we shall show that the super-ohmic case {s > 1) is bereft of any critical 
points. Finally in Sec. 15.31 we shall analyze the critical behavior when the bath 
spectrum is sub-ohmic in character (0 < s < 1). 



5.1 The RG equations for the ohmic case (s = 1) 

The situation of the ohmic bath spectrum is probably one of the most well under- 
stood case in the study of SBM. This is due to the fact that a linear dispersion of the 
bath DoS lends itself to an exact mapping to the anisotropic Kondo Model (AKM) 
13J |9). Due to this mapping, it is known that the critical dissipation occurs at a t = 1 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Fig. 6 RG flow for the ohmic SBM (s = 1). Here, the flow is constructed numerically by giving 
various initial (bare) values of the coupling constants /; and a. See discussion in Sec. l5.1l for the 
interpretation. 

for small non-zero A, and that the phase transition is of the Kosterlitz-Thouless type 
(infinite order). The flow equations that are found through the mapping to the AKM 
match the /3 functions that we obtained by using the Majorana representation (see 
Eq. [Pfl and Eq. fT3T >. From the structure of these j3 functions of the ohmic SBM, it 
is amply clear that the term —ah 2 drives the dissipative coupling to zero whenever 
a < 1. However, in the regime a > 1, it is now the transverse field term h that is 
driven to zero, with the dissipative coupling a renormalizing to a finite value. Fur- 
thermore, in the limit a > 1 we see that the RG equations, Eq.[l4]and Eq.[l5] have 
in fact a line of stable fixed points at zero field, the typical signature of a phase 
transition of the Kosterlitz-Thouless type. This discussion is encapsulated by Fig. [6] 
which represents the various RG trajectories that are obtained by numerically solv- 
ing Eq.[l4]and Eq.[l5] for various initial values of a and h. From this flow diagram, 
it is clear that there exists a separatrix such that for any value of a and /; that lies 
below the separatrix, the RG flow terminates at the line of fixed points, whereas if 
one were to start with initial value of a and h lying above the separatrix, the flow 
maintains the system in the de-localized phase. 
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Fig. 7 Flow equations of the super-ohmic SBM model (s > 1), starting with various initial values 
of a and h. For further details refer to Sec. l5.2l in the text. 

5.2 The RG equations for the super-ohmic case (s > I) 

In the situation where the bath spectrum is super-ohmic, i.e. s > 1, it can be readily 
argued that the system supports no critical fixed points. This fact can be essentially 
gleaned from solving the set of equations, Eq.Q~6]and Eq.[l5]numerically for various 
initial configurations of a and h, giving the results depicted in Fig. [7] One sees that 
for any initial value of the dissipation and the transverse field, the couplings always 
flow towards the limit h = °° and a = 0. This implies that for the super-ohmic case 
one always ends up in the de-localized phase, and no quantum phase trantion is 
allowed. 



5.3 The RG equations for the sub-ohmic case (0 < s < 1) 

Now, we turn our attention to the most interesting case, namely the one where the 
bath spectrum is endowed with a sub-ohmic dispersion, i.e. < s < 1. Since the 
mapping of the SBM to the AKM is invalidated in the sub-ohmic regime, the sit- 
uation of the SBM in the range < s < 1 was not fully appreciated until recent 
study from Numerical Renormalization Group (NRG) calculations J2], where a sec- 
ond order quantum phase transition was explicitely demonstrated for all < s < 1 . 
At this juncture, it should be noted that this localization/delocalization phase tran- 
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Fig. 8 Flow for the sub-ohmic SBM. Once again, the RG trajectories are plotted for various initial 
values of the transverse field /; and the dissipative coupling a. For further details refer Sec. l5.3l in 
the text. 

sition found in Ref. (2) is missed by the various other analytical treatments of the 
sub-ohmic SBM, such as variational ansatz or diagonalizations by unitary transfor- 
mations, but is correctly predicted by the flow equations derived above. The result- 
ing flow is plotted in Fig. [8] with a fixed point occuring at a c = 1 and /; = vl— s, 
perturbatively controlled for values of s close to 1 . 



6 Mapping to a long-ranged Ising Model 

In this section, we shall attempt to derive an effective model for the quantum phase 
transition discussed previously, based purely in terms of the bosonic mode <j>. This 
can be done by representing the spin now in terms of Abrikosov fermions |[T4ll . and 
then integrating the fermionic degrees of freedom perturbatively in A. By using this 
route we will see that the SBM model can be mapped to a 4 model with 0(1) 
symmetry and long-ranged interactions in imaginary time. The mapping works as 
follows: 

S=I/^/, (17) 

aa' 

Here, in Eq. [TT] the fl field are canonical fermions with an imaginary chemical 
potential 03), that redefines the Matsubara frequencies 0)„ — > co„ + nT /4, and a 
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are the three Pauli matrices. By using Eq. [17] to the defining Hamiltonian of the 
SBM, (Eq.0 can be re-written in terms of the following action: 



S = 



dxf 



/ + / dT^idr + Oi)^. (18) 



In the above equation / is a two component vector whose Hermitian conjugate is 
given by = (fl,fj). Also in Eq.[T8] the d t term is a consequence of time slicing 
when going into the path integral representation. The philosophy is now to formally 
integrate out the fermions to get an perturbative expansion in A of the effective 
action. This methodology of integrating out the fermions is very similar in spirit to 
the treatment by Hertz of the itinerant ferromagnet [7j, that we have already alluded 
to in the introductory section Q] This technical step can be formally performed as 
the fermionic sector is purely Gaussian, so that the effective theory reads: 



Seff — / dz^a' i (d T + (Qj)ai — Tr In 



A k 7X -,, 
-ff x +-(7-X(a 1 -- 



(19) 



Defining the "local mode" = (aj + fl,), the bath can be exactly encapsulated by 
the following Gaussian action, written with the Matsubara frequency v„ = 2nnT: 



Se G fT s = - o £%->„)#v n )<K-*v n ) 

P v„ 



The quantity % in the above equation is given by 

1 



V„)=L 



1 



jV n -COi -lV„-COj 

which can be re-expressed in terms of a spectral representation as: 

P(fl») 



(iVn) 



dco- 



Here, in Eq. 



iv„ - (0 

the bosonic density of states, p (©) is given by 



CO 



p(t») = — r— 6(co c + co)9((O c -(o)sgn(co). 



(20) 



(21) 



(22) 



(23) 



The Green's function % can be calculated by substituting the functional form of 
p(co) from Eq. [23] into Eq. [22] and then performing the integration over the fre- 
quency variable to. Once we have performed the integration, the resultant expres- 
sion can be easily inverted to obtain the functional form for 5f ~ 1 which is given in 
the low frequency limit by 
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2 v 2 / sin^eDc 

Thus, substituting the above functional form of the $# ~ into Eq. [20] we see that 
the Gaussian part of the action for the bath of harmonic oscillators, Eq. [19] can be 
written as: 



r-Gauss 1 Y" 1 

r v„ 



2 ' \ 2 ) sm¥(Q s + l 



|0(v„)| 2 (25) 



Now, that we have taken care of the Gaussian bath term in Eq. [19] we turn our 
attention to the 7>ln term, which can be Taylor expanded to obtain: 

°° 1 / A \ 2 " 

£-7>^-) , (26) 

n= 1 \ / 

wherein in the above equation, the quantity Go is endowed the functional form, 
Go = "°" l2 +T g ^ w here I2 is the usual 2x2 identity matrix. Now, we proceed to 

calculate the traces implicit in Eq.[26] The resultant expression, up to order A 4 , is 
then combined with Eq.[25j to obtain: 

Seff = J ^(r+A|v|')|*(iv)| 2 + J \/t«(*(t)) 4 . (27) 

Here r — s(O c /2 — A 2 / (4A) is a mass term that controls the distance to criticality, 
u = A 4 /(164 3 ) is the leading interaction term, and A « fi)* _1 . We note that this 
action is equivalent to an Ising model in imaginary time, with interaction decaying 
as 1/(t— %') , as expected from the quantum/classical equivalence 15] [2. Now, 
one can use simple power-counting arguments to capture the critical behavior of the 
long-ranged Ising model displayed in Eq.[27] thereby also understanding the critical 
behavior of the underlying microscopic model, Eq. [2] By doing a power counting 
analysis around the Gaussian fixed point one finds that the scale dimension of the 
O(0 4 ) term is [it] = 2s — 1. This implies that the for all s < 1/2 the scale dimension 
is negative thus implying that the critical behavior is mean field like. However, for 
s > 1/2 one needs to account for higher loop effects to capture the true critical 
behavior, leading to non trivial exponents with respect to the mean field values. 



7 A special identity 

In this section, we will derive a special identity that helps us to calculate the 
fully dressed bosonic propagator in terms of the spin-spin correlator &(t) = 
(<J : (t)o z (0)). We start with an effective action which is a variant of the one that 
can be obtained from Eq. [2] Thus, we have 
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+ £ f dTj a (T)a a (T)+ f dxdx'% l (x-x')<t>{x)<t>{x'). (28) 

a J J 

Here, in Eq. [28] the term SBen-y is the so-called Berry action that encodes the im- 
purity spin commutation relations. This term is not explicitly written down as its 
functional form relies on spin-coherent states, the discussion of which is beyond the 
scope of this manuscript. Also, in Eq. [28] J is a source term for the spin dynamics, 
and % is again the bare bosonic propagator. The spin-spin correlator can be easily 
derived from Eq. [28] by performing an appropriate functional differentiation of the 
partition function Z with respect to the source field J. Thus, we have 

^ t)= L/(t)T/(o) ^° (29) 

Now, in performing the technical calculations inherent in Eq.[29] one can re-express 
the bosonic field in terms of a new field = j- + 0. In doing so we use the fact 
that the partition function Z remains invariant under such a redefinition of the (f) 
field. Thus, re-expressing the partition function in terms of the fields and then 
performing the functional differentiation, we are led to the following relation that 
connects % z {x) to the full bosonic Green's function %(t— x') = (0(t)0(V)) and 
the bare bosonic Green's function 

Xz (T) = -±%-\T) + ± J dX l dX 2 % l {x l ~X)%\x 2 )^{x l )^(x 2 )) (30) 

By going into the frequency domain representation we can compactly re-write 
Eq.[30]as 

Xz^v n )^- TJ -j-—^ + T j-^- / — -j (31) 
X 2 %(iv n ) X 2 %(i Vn ) 2 

This identity couples the single particle full bosonic propagator Sfy to the spin-spin 
susceptibility % z , naively a four operators correlation function (see e.g. the decom- 
position onto Abrikosov fermions), and shows that both the bosonic field and the 
longitudinal spin density must become critical altogether at the quantum phase tran- 
sition. This formula becomes extremely useful in the context of diagrammatic ex- 
pansions that use the Majorana representation (introduced in section O, because 
the spin susceptibility is simply related to the single particle Majorana propagators, 
leading to a very powerful Ward identity. These further theoretical developments go 
however much beyond the scope of this review. 

However, in the next section, Sec.[H] we shall show the usefulness of the identity 
derived in this section to recover a well known result in the context of SBM models, 
the so-called Shiba's relation. 
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8 Shiba's Relation for the sub-ohmic spin boson model 



In this section we establish the Shiba's relation, usually discussed for the ohmic 
SBM, in the case of the sub-ohmic model s < 1 . This relation essentially connects 
the spin correlations at equilibrium to the zero-frequency spin susceptibility (via the 
free bath spectrum). To derive this we use the fact that the exact bosonic Green's 
function, on the real frequency axis, has the following low-frequency form (as can 
be checked by simple perturbative calculations from the effective action Eq. 127b 

S^(v) = (m + a°|v|' ' + i&°|v|'sgii(v)f 1 (32) 

where m is the renormalized mass driving the transition, and a s ,b s are non critical 
numerical coefficients. In the limit of small frequencies, the above reduces to: 

1 b° 

S^(v) = --i-^|v|'sgn(v) (33) 

Now, the identity Eq. |5T|gives us a relation that connects this full bosonic Green's 
function % to the spin susceptibility. At low frequency, and introducing the bare 
mass mo = l/%(0), one obviously gets: 



1 1 X 2 ,, x 

X'M (34) 

o 



m mo 4ml 



and 

x «( v ) = - b ^\ v f sg n(v)\x' z (0)} 2 (35) 

In Eq.|34]and Eq.[35] the quantities %\ an d Xz &K tne rea ^ an d imaginary part of the 
longitudinal spin susceptibility % z . The imaginary part of the bare bosonic Green's 

function reads: ^q'(v) = — ^jr^-sgn(v) = — ^j|v|' s sgn(v). Thus, substituting for 
in Eq.|35] we get: 

^'(v) = iy(|v|)sgn(v)[ Z ;(0)] 2 (36) 

Finally, from the definition that at T = the imaginary part of the spin suscep- 
tibility is related to the spin correlation function C(v) via the simple relation 
X" = sgn(v)C(v), we obtain 

C(v) = i/(|v|M0)] 2 (37) 



which is the generalized Shiba relation for the sub-ohmic spin boson model, and is 
valid in the low frequency limit for the whole delocalized phase. 



18 S. Florens, D. Venturelli and R. Narayanan 

9 On the Quantum to Classical mapping 



From the general arguments given in the introduction, and the detailed derivation of 
the classical effective theory Eq.|27]for the specific case of the spin boson model, 
the results of Ref. 11241 came as a surprise, since critical exponents associated to the 
spin magnetization (S z ) were numerically found by these authors to deviate from 
the expected mean field result for < s < 1/2. At the time of writing, this issue is 
still debated l25ll . see Ref. |6) for a more recent update on the question. 
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